Transcriptomic analysis of Anopheles gambiae from Benin reveals overexpression of salivary and cuticular proteins associated with cross-resistance to pyrethroids and organophosphates.

BACKGROUND: Insecticide resistance (IR) is one of the major threats to malaria vector control programs in endemic countries. However, the mechanisms underlying IR are poorly understood. Thus, investigating gene expression patterns related to IR can offer important insights into the molecular basis of IR in mosquitoes. In this study, RNA-Seq was used to characterize gene expression in Anopheles gambiae surviving exposure to pyrethroids (deltamethrin, alphacypermethrin) and an organophosphate (pirimiphos-methyl). RESULTS: Larvae of An. gambiae s.s. collected from Bassila and Djougou in Benin were reared to adulthood and phenotyped for IR using a modified CDC intensity bottle bioassay. The results showed that mosquitoes from Djougou were more resistant to pyrethroids (5X deltamethrin: 51.7% mortality; 2X alphacypermethrin: 47.4%) than Bassila (1X deltamethrin: 70.7%; 1X alphacypermethrin: 77.7%), while the latter were more resistant to pirimiphos-methyl (1.5X: 48.3% in Bassila and 1X: 21.5% in Djougou). RNA-seq was then conducted on resistant mosquitoes, non-exposed mosquitoes from the same locations and the laboratory-susceptible An. gambiae s.s. Kisumu strain. The results showed overexpression of detoxification genes, including cytochrome P450s (CYP12F2, CYP12F3, CYP4H15, CYP4H17, CYP6Z3, CYP9K1, CYP4G16, and CYP4D17), carboxylesterase genes (COEJHE5E, COE22933) and glutathione S-transferases (GSTE2 and GSTMS3) in all three resistant mosquito groups analyzed. Genes encoding cuticular proteins (CPR130, CPR10, CPR15, CPR16, CPR127, CPAP3-C, CPAP3-B, and CPR76) were also overexpressed in all the resistant groups, indicating their potential role in cross resistance in An. gambiae. Salivary gland protein genes related to 'salivary cysteine-rich peptide' and 'salivary secreted mucin 3' were also over-expressed and shared across all resistant groups. CONCLUSION: Our results suggest that in addition to metabolic enzymes, cuticular and salivary gland proteins could play an important role in cross-resistance to multiple classes of insecticides in Benin. These genes warrant further investigation to validate their functional role in An. gambiae resistance to insecticides.

CRISPR/Cas9-Mediated Knockout of the PxJHBP Gene Resulted in Increased Susceptibility to Bt Cry1Ac Protoxin and Reduced Lifespan and Spawning Rates in Plutella xylostella.

Juvenile hormone binding protein (JHBP) is a key regulator of JH signaling, and crosstalk between JH and 20-hydroxyecdysone (20E) can activate and fine-tune the mitogen-activated protein kinase cascade, leading to resistance to insecticidal proteins from Bacillis thuringiensis (Bt). However, the involvement of JHBP in the Bt Cry1Ac resistance of Plutella xylostella remains unclear. Here, we cloned a full-length cDNA encoding JHBP, and quantitative real-time PCR (qPCR) analysis showed that the expression of the PxJHBP gene in the midgut of the Cry1Ac-susceptible strain was significantly higher than that of the Cry1Ac-resistant strain. Furthermore, CRISPR/Cas9-mediated knockout of the PxJHBP gene significantly increased Cry1Ac susceptibility, resulting in a significantly shorter lifespan and reduced fertility. These results demonstrate that PxJHBP plays a critical role in the resistance to Cry1Ac protoxin and in the regulation of physiological metabolic processes associated with reproduction in adult females, providing valuable insights to improve management strategies of P. xylostella.

Aberrant splicing of a nicotinic acetylcholine receptor alpha 6 subunit is associated with spinosad tolerance in the thrips predator Orius laevigatus.

Susceptibility to insecticides is one of the limiting factors preventing wider adoption of natural enemies to control insect pest populations. Identification and selective breeding of insecticide tolerant strains of commercially used biological control agents (BCAs) is one of the approaches to overcome this constraint. Although a number of beneficial insects have been selected for increased tolerance to insecticides the molecular mechanisms underpinning these shifts in tolerance are not well characterised. Here we investigated the molecular mechanisms of enhanced tolerance of a lab selected strain of Orius laevigatus (Fieber) to the commonly used biopesticide spinosad. Transcriptomic analysis showed that spinosad tolerance is not a result of overexpressed detoxification genes. Molecular analysis of the target site for spinosyns, the nicotinic acetylcholine receptor (nAChR), revealed increased expression of truncated transcripts of the nAChR α6 subunit in the spinosad selected strain, a mechanism of resistance which was described previously in insect pest species. Collectively, our results demonstrate the mechanisms by which some beneficial biological control agents can evolve insecticide tolerance and will inform the development and deployment of insecticide-tolerant natural enemies in integrated pest management strategies.

Signatures of adaptation at key insecticide resistance loci in Anopheles gambiae in Southern Ghana revealed by reduced-coverage WGS.

Resistance to insecticides and adaptation to a diverse range of environments present challenges to Anopheles gambiae s.l. mosquito control efforts in sub-Saharan Africa. Whole-genome-sequencing is often employed for identifying the genomic basis underlying adaptation in Anopheles, but remains expensive for large-scale surveys. Reduced coverage whole-genome-sequencing can identify regions of the genome involved in adaptation at a lower cost, but is currently untested in Anopheles mosquitoes. Here, we use reduced coverage WGS to investigate population genetic structure and identify signatures of local adaptation in Anopheles mosquitoes across southern Ghana. In contrast to previous analyses, we find no structuring by ecoregion, with Anopheles coluzzii and Anopheles gambiae populations largely displaying the hallmarks of large, unstructured populations. However, we find signatures of selection at insecticide resistance loci that appear ubiquitous across ecoregions in An. coluzzii, and strongest in forest ecoregions in An. gambiae. Our study highlights resistance candidate genes in this region, and validates reduced coverage WGS, potentially to very low coverage levels, for population genomics and exploratory surveys for adaptation in Anopheles taxa.

Mutations and intron polymorphisms in voltage-gated sodium channel genes of different geographic populations of Culex pipiens pallens/Culex pipiens quinquefasciatus in China.

BACKGROUND: Culex pipiens pallens and Culex pipiens quinquefasciatus are the dominant species of Culex mosquitoes in China and important disease vectors. Long-term use of insecticides can cause mutations in the voltage-gated sodium channel (vgsc) gene of mosquitoes, but little is known about the current status and evolutionary origins of vgsc gene in different geographic populations. Therefore, this study aimed to determine the current status of vgsc genes in Cx. p. pallens and Cx. p. quinquefasciatus in China and to investigate the evolutionary inheritance of neighboring downstream introns of the vgsc gene to determine the impact of insecticides on long-term evolution. METHODS: Sampling was conducted from July to September 2021 in representative habitats of 22 provincial-level administrative divisions in China. Genomic DNA was extracted from 1308 mosquitoes, the IIS6 fragment of the vgsc gene on the nerve cell membrane was amplified using polymerase chain reaction, and the sequence was used to evaluate allele frequency and knockdown resistance (kdr) frequency. MEGA 11 was used to construct neighbor-joining (NJ) tree. PopART was used to build a TCS network. RESULTS: There were 6 alleles and 6 genotypes at the L1014 locus, which included the wild-type alleles TTA/L and CTA/L and the mutant alleles TTT/F, TTC/F, TCT/S and TCA/S. The geographic populations with a kdr frequency less than 20.00% were mainly concentrated in the regions north of 38° N, and the geographic populations with a kdr frequency greater than 80.00% were concentrated in the regions south of 30° N. kdr frequency increased with decreasing latitude. And within the same latitude, the frequency of kdr in large cities is relatively high. Mutations were correlated with the number of introns. The mutant allele TCA/S has only one intron, the mutant allele TTT/F has three introns, and the wild-type allele TTA/L has 17 introns. CONCLUSIONS: Cx. p. pallens and Cx. p. quinquefasciatus have developed resistance to insecticides in most regions of China. The neighboring downstream introns of the vgsc gene gradually decreased to one intron with the mutation of the vgsc gene. Mutations may originate from multiple mutation events rather than from a single origin, and populations lacking mutations may be genetically isolated.

A chromosome-level genome assembly of the soybean pod borer: insights into larval transcriptional response to transgenic soybean expressing the pesticidal Cry1Ac protein.

BACKGROUND: Genetically modified (GM) crop plants with transgenic expression of Bacillus thuringiensis (Bt) pesticidal proteins are used to manage feeding damage by pest insects. The durability of this technology is threatened by the selection for resistance in pest populations. The molecular mechanism(s) involved in insect physiological response or evolution of resistance to Bt is not fully understood. RESULTS: To investigate the response of a susceptible target insect to Bt, the soybean pod borer, Leguminivora glycinivorella (Lepidoptera: Tortricidae), was exposed to soybean, Glycine max, expressing Cry1Ac pesticidal protein or the non-transgenic parental cultivar. Assessment of larval changes in gene expression was facilitated by a third-generation sequenced and scaffolded chromosome-level assembly of the L. glycinivorella genome (657.4 Mb; 27 autosomes + Z chromosome), and subsequent structural annotation of 18,197 RefSeq gene models encoding 23,735 putative mRNA transcripts. Exposure of L. glycinivorella larvae to transgenic Cry1Ac G. max resulted in prediction of significant differential gene expression for 204 gene models (64 up- and 140 down-regulated) and differential splicing among isoforms for 10 genes compared to unexposed cohorts. Differentially expressed genes (DEGs) included putative peritrophic membrane constituents, orthologs of Bt receptor-encoding genes previously linked or associated with Bt resistance, and those involved in stress responses. Putative functional Gene Ontology (GO) annotations assigned to DEGs were significantly enriched for 36 categories at GO level 2, respectively. Most significantly enriched cellular component (CC), biological process (BP), and molecular function (MF) categories corresponded to vacuolar and microbody, transport and metabolic processes, and binding and reductase activities. The DEGs in enriched GO categories were biased for those that were down-regulated (≥ 0.783), with only MF categories GTPase and iron binding activities were bias for up-regulation genes. CONCLUSIONS: This study provides insights into pathways and processes involved larval response to Bt intoxication, which may inform future unbiased investigations into mechanisms of resistance that show no evidence of alteration in midgut receptors.

Whole transcriptomic analysis reveals overexpression of salivary gland and cuticular proteins genes in insecticide-resistant Anopheles arabiensis from Western Kenya.

BACKGROUND: Effective vector control is key to malaria prevention. However, this is now compromised by increased insecticide resistance due to continued reliance on insecticide-based control interventions. In Kenya, we have observed heterogenous resistance to pyrethroids and organophosphates in Anopheles arabiensis which is one of the most widespread malaria vectors in the country. We investigated the gene expression profiles of insecticide resistant An. arabiensis populations from Migori and Siaya counties in Western Kenya using RNA-Sequencing. Centers for Disease Control and Prevention (CDC) bottle assays were conducted using deltamethrin (DELTA), alphacypermethrin (ACYP) and pirimiphos-methyl (PMM) to determine the resistance status in both sites. RESULTS: Mosquitoes from Migori had average mortalities of 91%, 92% and 58% while those from Siaya had 85%, 86%, and 30% when exposed to DELTA, ACYP and PMM, respectively. RNA-Seq analysis was done on pools of mosquitoes which survived exposure ('resistant'), mosquitoes that were not exposed, and the insecticide-susceptible An. arabiensis Dongola strain. Gene expression profiles of resistant mosquitoes from both Migori and Siaya showed an overexpression mainly of salivary gland proteins belonging to both the short and long form D7 genes, and cuticular proteins (including CPR9, CPR10, CPR15, CPR16). Additionally, the overexpression of detoxification genes including cytochrome P450s (CYP9M1, CYP325H1, CYP4C27, CYP9L1 and CYP307A1), 2 carboxylesterases and a glutathione-S-transferase (GSTE4) were also shared between DELTA, ACYP, and PMM survivors, pointing to potential contribution to cross resistance to both pyrethroid and organophosphate insecticides. CONCLUSION: This study provides novel insights into the molecular basis of insecticide resistance in An. arabiensis in Western Kenya and suggests that salivary gland proteins and cuticular proteins are associated with resistance to multiple classes of insecticides.

Mesoporous silica nanospheres-mediated insecticide and antibiotics co-delivery system for synergizing insecticidal toxicity and reducing environmental risk of insecticide.

Mesoporous silica nanoparticles (MSNs) are efficient carriers of drugs, and are promising in developing novel pesticide formulations. The cotton aphids Aphis gossypii Glover is a world devastating insect pest. It has evolved high level resistance to various insecticides thus resulted in the application of higher doses of insecticides, which raised environmental risk. In this study, the MSNs based pesticide/antibiotic delivery system was constructed for co-delivery of ampicillin (Amp) and imidacloprid (IMI). The IMI@Amp@MSNs complexes have improved toxicity against cotton aphids, and reduced acute toxicity to zebrafish. From the 16S rDNA sequencing results, Amp@MSNs, prepared by loading ampicillin to the mesoporous of MSNs, greatly disturbed the gut community of cotton aphids. Then, the relative expression of at least 25 cytochrome P450 genes of A. gossypii was significantly suppressed, including CYP6CY19 and CYP6CY22, which were found to be associated with imidacloprid resistance by RNAi. The bioassay results indicated that the synergy ratio of ampicillin to imidacloprid was 1.6, while Amp@MSNs improved the toxicity of imidacloprid by 2.4-fold. In addition, IMI@Amp@MSNs significantly improved the penetration of imidacloprid, and contributed to the amount of imidacloprid delivered to A. gossypii increased 1.4-fold. Thus, through inhibiting the relative expression of cytochrome P450 genes and improving penetration of imidacloprid, the toxicity of IMI@Amp@MSNs was 6.0-fold higher than that of imidacloprid. The greenhouse experiments further demonstrated the enhanced insecticidal activity of IMI@Amp@MSNs to A. gossypii. Meanwhile, the LC50 of IMI@Amp@MSNs to zebrafish was 3.9-fold higher than that of IMI, and the EC50 for malformation was 2.8-fold higher than IMI, respectively, which indicated that the IMI@Amp@MSNs complexes significantly reduced the environmental risk of imidacloprid. These findings encouraged the development of pesticide/antibiotic co-delivery nanoparticles, which would benefit pesticide reduction and environmental safety.

The frequency of kdr and ace-1 alleles in Anopheles gambiae s.l. before and during indoor residual spraying (IRS) implementation and four years after IRS withdrawal in three districts in Atacora, Benin.

BACKGROUND: Indoor residual spraying (IRS) was first implemented in the Atacora department, Benin from 2011 to 2012 using bendiocarb (carbamate) followed by annual spraying with pirimiphos-methyl (organophosphate) from 2013 to 2018. Before and after IRS implementation in Atacora, standard pyrethroid insecticide-treated bed nets were the main method of vector control in the area. This study investigated the knockdown resistance (kdr) gene (L1014F) and the acetylcholinesterase (ace-1) gene (G119S), before and during IRS implementation, and 4-years after IRS withdrawal from Atacora. This was done to assess how changes in insecticide pressure from indoor residual spraying may have altered the genotypic resistance profile of Anopheles gambiae s.l. METHOD: Identification of sibling species of An. gambiae s.l. and detection of the L1014F mutation in the kdr gene and G119S mutation in ace-1 genes was done using molecular analysis. Allelic and genotypic frequencies were calculated and compared with each other before and during IRS implementation and 4 years after IRS withdrawal. The Hardy-Weinberg equilibrium and genetic differentiation within and between populations were assessed. RESULTS: Prevalence of the L1014F mutation in all geographic An. gambiae s.l. (An. gambiae s.s., Anopheles. coluzzii, Anopheles. arabiensis, and hybrids of "An. gambiae s.s. and An. coluzzii") populations increased from 69% before IRS to 87% and 90% during and after IRS. The G119S allele frequency during IRS (20%) was significantly higher than before IRS implementation (2%). Four years after IRS withdrawal, allele frequencies returned to similar levels as before IRS (3%). Four years after IRS withdrawal, the populations showed excess heterozygosity at the ace-1 gene and deficit heterozygosity at the kdr gene, whereas both genes had excess heterozygosity before and during IRS (FIS < 0). No genetic differentiation was observed within the populations. CONCLUSIONS: This study shows that the withdrawal of IRS with bendiocarb and pirimiphos-methyl may have slowed down the selection of individual mosquitoes with ace-1 resistance alleles in contrast to populations of An. gambiae s.l. with the L1014F resistance allele of the kdr gene. This may suggest that withdrawing the use of carbamates or organophosphates from IRS or rotating alternative insecticides with different modes of action may slow the development of ace-1 insecticide-resistance mutations. The increase in the prevalence of the L1014F mutation of the kdr gene in the population, despite the cessation of IRS, could be explained by the growing use of pyrethroids and DDT in agriculture and for other domestic use. More observational studies in countries where carbamates or organophosphates are still being used as public health insecticides may provide additional insights into these associations.

The impact of agrochemical pollutant mixtures on the selection of insecticide resistance in the malaria vector Anopheles gambiae: insights from experimental evolution and transcriptomics.

BACKGROUND: There are several indications that pesticides used in agriculture contribute to the emergence and spread of resistance of mosquitoes to vector control insecticides. However, the impact of such an indirect selection pressure has rarely been quantified and the molecular mechanisms involved are still poorly characterized. In this context, experimental selection with different agrochemical mixtures was conducted in Anopheles gambiae. The multi-generational impact of agrochemicals on insecticide resistance was evaluated by phenotypic and molecular approaches. METHODS: Mosquito larvae were selected for 30 generations with three different agrochemical mixtures containing (i) insecticides, (ii) non-insecticides compounds, and (iii) both insecticide and non-insecticide compounds. Every five generations, the resistance of adults to deltamethrin and bendiocarb was monitored using bioassays. The frequencies of the kdr (L995F) and ace1 (G119S) target-site mutations were monitored every 10 generations. RNAseq was performed on all lines at generation 30 in order to identify gene transcription level variations and polymorphisms associated with each selection regime. RESULTS: Larval selection with agrochemical mixtures did not affect bendiocarb resistance and did not select for ace1 mutation. Contrastingly, an increased deltamethrin resistance was observed in the three selected lines. Such increased resistance was not majorly associated with the presence of kdr L995F mutation in selected lines. RNA-seq identified 63 candidate resistance genes over-transcribed in at least one selected line. These include genes coding for detoxification enzymes or cuticular proteins previously associated with insecticide resistance, and other genes potentially associated with chemical stress response. Combining an allele frequency filtering with a Bayesian FST-based genome scan allowed to identify genes under selection across multiple genomic loci, supporting a multigenic adaptive response to agrochemical mixtures. CONCLUSION: This study supports the role of agrochemical contaminants as a significant larval selection pressure favouring insecticide resistance in malaria vectors. Such selection pressures likely impact kdr mutations and detoxification enzymes, but also more generalist mechanisms such as cuticle resistance, which could potentially lead to cross-tolerance to unrelated insecticide compounds. Such indirect effect of global landscape pollution on mosquito resistance to public health insecticides deserves further attention since it can affect the nature and dynamics of resistance alleles circulating in malaria vectors and impact the efficacy of control vector strategies.

Cross-pollination in seed-blended refuge and selection for Vip3A resistance in a lepidopteran pest as detected by genomic monitoring.

The evolution of pest resistance to management tools reduces productivity and results in economic losses in agricultural systems. To slow its emergence and spread, monitoring and prevention practices are implemented in resistance management programs. Recent work suggests that genomic approaches can identify signs of emerging resistance to aid in resistance management. Here, we empirically examined the sensitivity of genomic monitoring for resistance management in transgenic Bt crops, a globally important agricultural innovation. Whole genome resequencing of wild North American Helicoverpa zea collected from non-expressing refuge and plants expressing Cry1Ab confirmed that resistance-associated signatures of selection were detectable after a single generation of exposure. Upon demonstrating its sensitivity, we applied genomic monitoring to wild H. zea that survived Vip3A exposure resulting from cross-pollination of refuge plants in seed-blended plots. Refuge seed interplanted with transgenic seed exposed H. zea to sublethal doses of Vip3A protein in corn ears and was associated with allele frequency divergence across the genome. Some of the greatest allele frequency divergence occurred in genomic regions adjacent to a previously described candidate gene for Vip3A resistance. Our work highlights the power of genomic monitoring to sensitively detect heritable changes associated with field exposure to Bt toxins and suggests that seed-blended refuge will likely hasten the evolution of resistance to Vip3A in lepidopteran pests.

Monitoring Aedes populations for arboviruses, Wolbachia, insecticide resistance and its mechanisms in various agroecosystems in Benin.

Aedes mosquitoes are the main vectors of arboviruses in Benin. Cases of dengue have been reported in Benin with all four serotypes of the virus actively circulating in this region. Some agricultural settings are known to harbor Aedes vectors responsible for the transmission of arboviruses. The massive use of certain insecticides in agricultural settings has probably contributed to insecticide resistance in these vectors. In Benin, the susceptibility of arbovirus vectors to insecticides is poorly studied. In addition, the distribution of Wolbachia spp., which is used against some arboviruses is unknown. Moreover, there is limited information regarding the vectors responsible for the transmission of arboviruses in Benin. This present study monitored the species composition, arboviruses, and Wolbachia symbiont status, as well as the phenotypic and molecular insecticide resistance profile of Aedes populations from three agroecosystems in Benin. Aedes species identification was performed morphologically and confirmed using qPCR. (RT)-qPCR assay was applied for monitoring the presence of DENV, CHIKV, ZIKV, and WNV pathogens as well as for naturally occurring Wolbachia symbionts. Insecticide resistance was assessed phenotypically, by permethrin (0.75%) exposure of Adults (F0) using World Health Organization (WHO) bioassay protocols, and at the molecular level, using TaqMan (RT)-qPCR assays for assessing knock-down resistance (kdr) mutations (F1534C, V1016G/I, and S989P) and the expression levels of eight detoxification genes (P450s from the CYP9 and CYP6 families, carboxylesterases and glutathione-S-transferases). Aedes aegypti (Ae. aegypti) mosquitoes were the most abundant (93.9%) in the three agroecosystems studied, followed by Aedes albopictus (Ae. albopictus) mosquitoes (6.1%). No arboviruses were detected in the study's mosquito populations. Naturally occurring Wolbachia symbionts were present in 7 pools out of 15 pools tested. This could influence the effectiveness of vector control strategies based on exogenously introduced Wolbachia, all present in the three agroecosystems. Full susceptibility to permethrin was observed in all tested populations of Ae. albopictus. On the contrary, Ae. aegypti were found to be resistant in all three agroecosystem sites except for banana plantation sites, where full susceptibility was observed. Molecular analysis revealed that individual target site resistance kdr mutations F1534C and V1016G/I were detected in most Ae. aegypti populations. Additionally, double mutant (F1534C + V1016G/I) mosquitoes were found in some populations, and in one case, triple mutant (F1534C + V1016G/I + S989P) mosquitoes were detected. Metabolic resistance, as reflected by overexpression of three P450 genes (CYP6BB2, CYP9J26, and CYP9J32), was also detected in Ae. aegypti mosquitoes. Our study provides information that could be used to strategize future vector control strategies and highlights the importance of continuing vector surveillance. Future studies should assess the effect of piperonyl butoxide (PBO) on metabolic resistance and identify the different strains of Wolbachia spp., to choose the best vector control strategies in Benin.

Key role of the ryanodine receptor I4790K mutation in mediating diamide resistance in Plutella xylostella.

The diamondback moth Plutella xylostella, a global insect pest of cruciferous vegetables, has evolved resistance to many classes of insecticides including diamides. Three point mutations (I4790M, I4790K, and G4946E) in the ryanodine receptor of P. xylostella (PxRyR) have been identified to associate with varying levels of resistance. In this study, we generated a knockin strain (I4790K-KI) of P. xylostella, using CRISPR/Cas9 to introduce the I4790K mutation into PxRyR of the susceptible IPP-S strain. Compared to IPP-S, the edited I4790K-KI strain exhibited high levels of resistance to both anthranilic diamides (chlorantraniliprole 1857-fold, cyantraniliprole 1433-fold) and the phthalic acid diamide flubendiamide (>2272-fold). Resistance to chlorantraniliprole in the I4790K-KI strain was inherited in an autosomal and recessive mode, and genetically linked with the I4790K knockin mutation. Computational modeling suggests the I4790K mutation reduces the binding of diamides to PxRyR by disrupting key hydrogen bonding interactions within the binding cavity. The approximate frequencies of the 4790M, 4790K, and 4946E alleles were assessed in ten geographical field populations of P. xylostella collected in China in 2021. The levels of chlorantraniliprole resistance (2.3- to 1444-fold) in these populations were significantly correlated with the frequencies (0.017-0.917) of the 4790K allele, but not with either 4790M (0-0.183) or 4946E (0.017-0.450) alleles. This demonstrates that the PxRyR I4790K mutation is currently the major contributing factor to chlorantraniliprole resistance in P. xylostella field populations within China. Our findings provide in vivo functional evidence for the causality of the I4790K mutation in PxRyR with high levels of diamide resistance in P. xylostella, and suggest that tracking the frequency of the I4790K allele is crucial for optimizing the monitoring and management of diamide resistance in this crop pest.

Cytochrome P450 CYP6EM1 Underpins Dinotefuran Resistance in the Whitefly Bemisia tabaci.

Being a destructive pest worldwide, the whitefly Bemisia tabaci has evolved resistance to neonicotinoid insecticides. The third-generation neonicotinoid dinotefuran has commonly been applied to the control of the whitefly, but its underlying mechanism is currently unknown. On the base of our transcriptome data, here we aim to investigate whether the cytochrome P450 CYP6EM1 underlies dinotefuran resistance in the whitefly. Compared to the susceptible strain, the CYP6EM1 gene was found to be highly expressed in both laboratory and field dinotefuran-resistant populations. Upon exposure to dinotefuran, the mRNA levels of CYP6EM1 were increased. These results demonstrate the involvement of this gene in dinotefuran resistance. Loss and gain of functional studies in vivo were conducted through RNAi and transgenic Drosophila melanogaster assays, confirming the role of CYP6EM1 in conferring such resistance. In a metabolism assay in vitro, the CYP6EM1 protein could metabolize 28.11% of dinotefuran with a possible dinotefuran-dm-NNO metabolite via UPLC-QTOF/MS. Docking of dinotefuran to the CYP6EM1 protein showed a good binding affinity, with an energy of less than -6.0 kcal/mol. Overall, these results provide compelling evidence that CYP6EM1 plays a crucial role in the metabolic resistance of B. tabaci to dinotefuran. Our work provides new insights into the mechanism underlying neonicotinoid resistance and applied knowledge that can contribute to sustainable control of a global pest such as whitefly.

Mitochondrial coding genes mediate insecticide tolerance in the oriental fruit fly, Bactrocera dorsalis (Hendel).

The oriental fruit fly, Bactrocera dorsalis (Hendel), an invasive insect pest infesting fruits and vegetables, possesses a remarkable capacity for environmental adaptation. The investigation of behind mechanisms of the stress adaptability in B. dorsalis holds significantly practical relevance. Previous studies on the molecular mechanism underlying stress resistance in B. dorsalis have predominantly focused on nuclear-coding genes, with limited exploration on organelle-coding genes. In this study, we assessed alterations in the mitochondrial physiological parameters of B. dorsalis under exposure to malathion, avermectin, and beta-cypermethrin at LD50 dosages. The results showed that all three insecticides were capable of reducing mitochondrial complex IV activity and ATP content. Expression patterns of mitochondrial coding genes across different developmental stages, tissues and insecticide exposures were analyzed by RT-qPCR. The results revealed that these mitochondrial coding genes were expressed in various tissues and at different developmental stages. Particularly noteworthy, atp6, cox2, and cytb exhibited substantial up-regulation in response to malathion and avermectin treatment. Furthermore, RNAi-mediated knockdown of atp6 and cox2 resulted in the increased toxicity of malathion and avermectin against B. dorsalis, and cox2 silencing was also associated with the decreased complex IV activity. These findings suggest that atp6 and cox2 most likely play pivotal roles in mediating tolerance or resistance to malathion and avermectin in B. dorsalis. Our results provide novel insights into the role of mitochondrial coding genes in conferring tolerance to insecticides in B. dorsalis, with practical implications for controlling this pest in the field.

Individual transmembrane domains of SfABCC2 from Spodoptera frugiperda do not serve as functional Cry1F receptors.

The fall armyworm (Spodoptera frugiperda) is a major global pest causing severe damage to various crops, especially corn. Transgenic corn producing the Cry1F pesticidal protein from the bacterium Bacillus thuringiensis (Cry1F corn) showed effectiveness in controlling this pest until S. frugiperda populations at locations in North and South America evolved practical resistance. The mechanism for practical resistance involved disruptive mutations in an ATP binding cassette transporter subfamily C2 gene (SfABCC2), which serves as a functional Cry1F receptor in the midgut cells of susceptible S. frugiperda. The SfABCC2 protein contains two transmembrane domains (TMD1 and TMD2), each with a cytosolic nucleotide (ATP) binding domain (NBD1 and NBD2, respectively). Previous reports have demonstrated that disruptive mutations in TMD2 were linked with resistance to Cry1F, yet whether the complete SfABCC2 structure is needed for receptor functionality or if a single TMD-NBD protein can serve as functional Cry1F receptor remains unknown. In the present study, we separately expressed TMD1 and TMD2 with their corresponding NBDs in cultured insect cells and tested their Cry1F receptor functionality. Our results show that the complete SfABCC2 structure is required for Cry1F receptor functionality. Moreover, binding competition assays revealed that Cry1F specifically bound to SfABCC2, whereas neither SfTMD1-NBD1 nor SfTMD2-NBD2 exhibited any significant binding. These results provide insights into the molecular mechanism of Cry1F recognition by SfABCC2 in S. frugiperda, which could facilitate the development of more effective insecticidal proteins.

Knock-in mutagenesis in Drosophila Rdl underscores the critical role of the conserved M3 glycine in mediating the actions of broflanilide and isocycloseram on GABA receptors.

γ-Aminobutyric acid receptors (GABARs) are crucial targets for pest control chemicals, including meta-diamide and isoxazoline insecticides, which act as negative allosteric modulators of insect GABARs. Previous cell-based assays have indicated that amino acid residues in the transmembrane cavity between adjacent subunits of Drosophila RDL GABAR (i.e., Ile276, Leu280, and Gly335) are involved in mediating the action of meta-diamides. In this study, to confirm this result at the organismal level, we employed CRISPR/Cas9-mediated genome editing, generated six transgenic Drosophila strains carrying substitutions in these amino acid residues, and investigated their sensitivity to broflanilide and isocycloseram. Flies homozygous for the I276F mutation did not exhibit any change in sensitivity to the tested insecticides compared to the control flies. Conversely, I276C homozygosity was lethal, and heterozygous flies exhibited ∼2-fold lower sensitivity to broflanilide than the control flies. Flies homozygous for the L280C mutation survived into adulthood but exhibited infertility. Both heterozygous and homozygous L280C flies exhibited ∼3- and âˆ¼20-fold lower sensitivities to broflanilide and isocycloseram, respectively, than the control flies. The reduction in sensitivity to isocycloseram in L280C flies diminished to ∼3-fold when treated with piperonyl butoxide. Flies homozygous for the G335A mutation reached the adult stage. However, they were sterile, had small bodies, and exhibited reduced locomotion, indicating the critical role of Gly335 in RDL function. These flies exhibited markedly increased tolerance to topically applied broflanilide and isocycloseram, demonstrating that the conserved Gly335 is the target of the insecticidal actions of broflanilide and isocycloseram. Considering the significant fitness costs, the Gly335 mutation may not pose a serious risk for the development of resistance in field populations of insect pests. However, more careful studies using insect pests are needed to investigate whether our perspective applies to resistance development under field conditions.

Pyrethroid resistance status and co-occurrence of V1016G, F1534C and S989P mutations in the Aedes aegypti population from two dengue outbreak counties along the China-Myanmar border.

BACKGROUND: Over the past two decades, dengue fever (DF) has emerged as a significant arboviral disease in Yunnan province, China, particularly in the China-Myanmar border area. Aedes aegypti, an invasive mosquito species, plays a crucial role in transmitting the dengue virus to the local population. Insecticide-based vector control has been the primary tool employed to combat DF, but the current susceptibility status of Ae. aegypti to commonly used insecticides is unknown. Assessment of Ae. aegypti resistance to pyrethroid insecticides and an understanding of the underlying mechanisms of this resistance in the China-Myanmar border region is of significant strategic importance for effectively controlling the DF epidemic in the area. METHODS: Aedes aegypti larvae collected from Ruili and Gengma counties in Yunnan Province were reared to adults in the laboratory and tested for susceptibility to three pyrethroid insecticides (3.20% permethrin, 0.08% lambda-cyhalothrin and 0.20% deltamethrin) by the standard WHO susceptibility bioassay. Genotyping of mutations in the knockdown gene (kdr), namely S989P, V1016G and F1534C, that are responsible for resistance to pyrethroid insecticides was performed using allele-specific PCR methods. A possible association between the observed resistant phenotype and mutations in the voltage-gated sodium channel gene (VGSC) was also studied. RESULTS: Aedes aegypti mosquitoes collected from the two counties and reared in the laboratory were resistant to all of the pyrethroids tested, with the exception of Ae. aegypti from Gengma County, which showed sensitivity to 0.20% deltamethrin. The mortality rate of Ae. aegypti from Ruili county exposed to 3.20% permethrin did not differ significantly from that of Ae. aegypti from Gengma County (χ2 = 0.311, P = 0.577). By contrast, the mortality rate of Ae. aegypti from Ruili County exposed to 0.08% lambda-cyhalothrin and 0.20% deltamethrin, respectively, was significantly different from that of Ae. aegypti from Gengma. There was no significant difference in the observed KDT50 of Ae. aegypti from the two counties to various insecticides. Four mutation types and 12 genotypes were detected at three kdr mutation sites. Based on results from all tested Ae. aegypti, the V1016G mutation was the most prevalent kdr mutation (100% prevalence), followed by the S989P mutation (81.6%) and the F1534C mutation (78.9%). The constituent ratio of VGSC gene mutation types was significantly different in Ae. aegypti mosquitoes from Ruili and those Gengma. The triple mutant S989P + V1016G + F1534C was observed in 274 Ae. aegypti mosquitoes (60.8%), with the most common genotype being SP + GG + FC (31.4%). The prevalence of the F1534C mutation was significantly higher in resistant Ae. aegypti from Ruili (odds ratio [OR] 7.43; 95% confidence interval [CI] 1.71-32.29; P = 0.01) and Gengma (OR 9.29; 95% CI 3.38-25.50; P = 0.00) counties than in susceptible Ae. aegypti when exposed to 3.20% permethrin and 0.08% lambda-cyhalothrin, respectively. No significant association was observed in the triple mutation genotypes with the Ae. aegypti population exposed to 3.20% permethrin and 0.20% deltamethrin resistance (P > 0.05), except for Ae. aegypti from Gengma County when exposed to 0.08% lambda-cyhalothrin (OR 2.86; 95% CI 1.20-6.81; P = 0.02). CONCLUSIONS: Aedes aegypti from Ruili and Gengma counties have developed resistance to various pyrethroid insecticides. The occurrence of multiple mutant sites in VGSC strongly correlated with the high levels of resistance to pyrethroids in the Ae. aegypti populations, highlighting the need for alternative strategies to manage the spread of resistance. A region-specific control strategy for dengue vectors needs to be implemented in the future based on the status of insecticide resistance and kdr mutations.

Characterization of the inheritance of field-evolved resistance to diamides in the fall armyworm (Spodoptera frugiperda) (Lepidoptera: Noctuidae) population from Puerto Rico.

The fall armyworm (Spodoptera frugiperda) is one of the most destructive pests of corn. New infestations have been reported in the East Hemisphere, reaching India, China, Malaysia, and Australia, causing severe destruction to corn and other crops. In Puerto Rico, practical resistance to different mode of action compounds has been reported in cornfields. In this study, we characterized the inheritance of resistance to chlorantraniliprole and flubendiamide and identified the possible cross-resistance to cyantraniliprole and cyclaniliprole. The Puerto Rican (PR) strain showed high levels of resistance to flubendiamide (RR50 = 2,762-fold) and chlorantraniliprole (RR50 = 96-fold). The inheritance of resistance showed an autosomal inheritance for chlorantraniliprole and an X-linked inheritance for flubendiamide. The trend of the dominance of resistance demonstrated an incompletely recessive trait for H1 (♂ SUS × â™€ PR) × and an incompletely dominant trait for H2 (♀ SUS × â™‚ PR) × for flubendiamide and chlorantraniliprole. The PR strain showed no significant presence of detoxification enzymes (using synergists: PBO, DEF, DEM, and VER) to chlorantraniliprole; however, for flubendiamide the SR = 2.7 (DEM), SR = 3.2 (DEF) and SR = 7.6 (VER) indicated the role of esterases, glutathione S- transferases and ABC transporters in the metabolism of flubendiamide. The PR strain showed high and low cross-resistance to cyantraniliprole (74-fold) and cyclaniliprole (11-fold), respectively. Incomplete recessiveness might lead to the survival of heterozygous individuals when the decay of diamide residue occurs in plant tissues. These results highlight the importance of adopting diverse pest management strategies, including insecticide rotating to manage FAW populations in Puerto Rico and other continents.